Optomechanical cooling of levitated spheres with doubly-resonant fields 
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Optomechanical cooling of levitated dielectric particles represents a promising new approach in 
the quest to cool small mechanical resonators towards their quantum ground state. We investigate 
two-mode cooling of levitated nanospheres in a self-trapping regime. We identify a rich structure 
of split sidebands (by a mechanism unrelated to usual strong-coupling effects) and strong cooling 
even when one mode is blue detuned. We show the best regimes occur when both optical fields 
cooperatively cool and trap the nanosphere, where cooling rates are over an order of magnitude 
faster compared to corresponding single-sideband cooling rates. 

PACS numbers: 



Extraordinary progress has been made in the last half- 
dozen years [H, 0| towards the final goal of cooling a 
small mechanical resonator down to its quantum ground 
state and hence to realise quantum behavior in a macro- 
scopic system. Implementations include cavity cool- 
ing of micromirrors on cantilevers [3|-|6( , dielectric mem- 
branes in Fabry Perot cavities [7|; radial and whisper- 
ing gallery modes of optical microcavities [8] and nano- 
electromechanical systems [9(. Indeed the realizations 
span 12 orders of magnitude [2J, up to and including 
the LIGO gravity wave experiments. Corresponding ad- 
vances in the theory of optomechanical cooling have also 
been made |10l-fl3j. 

Over the last year or so, a promising new paradigm 
has been attracting much interest: several groups [14T - tl7j 
have now proposed schemes for optomechanical cooling of 
levitated dielectric particles, including nanospheres and 
The important advantage is the elimi- 



14] 



even viruses 

nation of the mechanical support, a dominant source of 

heating noise. In general, these proposals involve two 

fields, one for trapping and one for cooling. This may 

involve an optical cavity mode plus a separate trap; or 

two optical cavity modes, the so-called "self-trapping" 

scenario. 

Mechanical oscillators in the self-trapping regime differ 
from other optomechanically-cooled devices in a second 
fundamental respect (in addition to the absence of me- 
chanical support): the mechanical frequency, wjw, asso- 
ciated with centre of mass oscillations is not an intrinsic 
feature of the resonator but is determined by the opti- 
cal field. In particular, it is a function of one or both of 
the detuning frequencies, d± and 5%, of the optical modes. 
Cooling, in general, occurs when wm is resonantly red de- 
tuned with either of the detuning frequencies (i.e. neg- 
ative #1,2 is associated with cooling). For self-trapping 
systems, this means u>m(^Ii^2) ~ — <5i,2 so the relevant 
frequencies are not independent. 

The full implications of this nonlinear interdependence 
of the resonant frequencies have not yet been fully elu- 
cidated. We show here for the first time that it leads 
to a rich landscape of split sidebands. The mechanism 
here is unrelated to splittings seen in experiments in the 



strong-coupling regime |l9[. However, it results in ex- 
tremely favourable cooling regimes, where two (or more) 
cooling sidebands approach each other. We term this 
the "double-resonance" regime. We find it can produce 
cooling rates nearly two orders of magnitude stronger 
than the corresponding "single-resonance" case. A self- 




FIG. 1: Schematic set-up: a levitated nanosphere is trapped 
and cooled cooperatively by two optical modes. The optical 
potential for each mode is shown. 

trapping Hamiltonian was investigated in [lfj and corre- 
sponds to the set-up illustrated in Figfl] 

P 2 



t„ A ;,t„ : ' /l 5 t a , 

Aa\ai cos 2 k\X + Ei(a\ + Si) + REi(a\ + a 2 )(l) 
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Two optical field modes a\ y 2 are coupled to a nanosphere 
with centre of mass position x. H is given in the ro- 
tating frame of the laser which drives the modes with 
amplitudes E\ and RE\ respectively. In [15|], the phase 
between the optical potentials was chosen to be <\> = ir/4; 
the study focussed primarily on the Si ~ regime, where 
the ai mode is responsible exclusively for trapping while 
the 6,2 mode alone provides cooling. Previous studies 
[lj, [la, LL2I ah analysed mechanical oscillations about an 
equilibrium position x ~ 0, corresponding to the antin- 
ode of the trapping mode (field 1). Below, this scenario 
is referred to as the "single-resonance" regime. 

Here we investigate the effects of relaxing all these 
restrictions and find interesting and unexpected impli- 
cations. We take <j> = 7r/4; the cooling field is driven 



more weakly than the trapping field, but with a ratio 
R ~ 0.1 — 1. Below, our analytical expressions cover arbi- 
trary k, R, Ei, A, but we compare with an illustrative set 
of experimentally plausible parameters: we take a cavity 
damping k = 6x 10°Hz. We considered driving powers 

in the range P ~ 1 — 10 mW, where P = — 1 . For a 
laser of wavelength A = 1064 nm and a cavity of length 
L ~ 1 cm, waist 25/xm we consider a silica nanosphere of 
lOOnm radius and hence a coupling strength A~3x 10 5 
Hz. To obtain a dephasing of 7r/4 between the two modes 
near the centre of the cavity, the frequency difference be- 
tween the modes is |wi — W2I ~ 27r x 10 GHz. This far 
exceeds the detunings (5 12 ~ 1 MHz and also the mechan- 
ical frequencies %. Thus the photons are completely 
distinguishable and can be read out and driven sepa- 
rately. Nevertheless, since oj 12 ~ 10 14 Hz, we approx- 
imate fci ~ fc 2 = k. However, this situation is distinct 
from the ring-cavity proposal of [17j | where R = and 
mode 2 is undriven but is populated exclusively by scat- 
tering from mode 1; there, the photons are of precisely 
the same frequency and thus there is a single detuning 
parameter involved. 



Fig(2]illustrates the behavior for R = 0.5. It shows that 
allowing both fields to cooperatively trap and cool yields 
more than an additive improvement. We denote by rl 
and rl the set of detunings corresponding to cooling res- 
onances of modes 1 and 2 respectively: FigJ5] shows that 
these resonances unexpectedly split and separate into 
new cooling resonances rl± and r2±. These can over- 
lap, to give very strong cooling associated with multiple 
resonances. Below, we also show that the usually studied 
single-field resonance regime r2 can attain only a maxi- 
mal cooling rate T ~ ^j^; and that for strong driving, 
r oc \jE\\ thus, the cooling falls with increasing driving 
and it is hard to achieve optimal cooling regimes where 

r ~ A ~ k. For the double-field resonances, in contrast, 

1/2 
r oc Ei , the cooling increases with Ei and can more 

easily reach optimal cooling. Very strong cooling is ap- 
parent even for regimes where one mode is blue detuned. 
In addition, although there is no direct coupling between 
optical modes, double-resonances offer the prospect of 
strong (albeit second-order) coupling and entangling of 
the two modes via the nanosphere, within a single cav- 
ity. This includes simultaneous resonant/antiresonant 
regimes (indicated by the crossing of r2+ and al— in 
FigJ2|) where one mode resonantly heats, while the other 
resonantly cools the mechanical mode. 



The dynamics depends on k, m, A, k, Si, S2,Ei,R. 
However transforming to scaled variables reduces this 
complexity. We rescale position, time and field variables 
as follows: kx —tx,At—t i, then ai, 2 — > ^¥&\,i Note 
that below we drop all the tildes but it is implicit that 
all variables are scaled in the resulting Heisenberg equa- 
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FIG. 2: (Colour online)Upper Panels: comparison between 
numerical optical cooling rates (without linearisation) and an 
analytical expression (Eq[5} from linearised dynamics, show- 
ing excellent agreement. R = 0.5. The curves corresponding 
to different values of 5i ± 2MHz are shifted relative to each 
other. At single-resonance r2, field 2 is resonant with the 
oscillator and is exclusively responsible for cooling; field 1 is 
resonant with the cavity (Af = 0) and traps the sphere. Sub- 
sequently r2± appear: they are cooling resonances of field 
2, split by field 1; conversely rl± are cooling resonances of 
field 1, split by field 2. r2— , rl— , rl+ overlap giving a broad 
region of very strong cooling. The al± are heating (Stokes) 
resonances of field 1. The al— can coincide with the cool- 
ing resonance (r2+). Here field 2 absorbs phonons as fast 
as field 1 emits them. Lower Panels: Show the unshifted 
cooling curves at the single-field cooling resonance r2 and 
double- resonant cooling, showing that the latter gives over an 
order of magnitude stronger cooling. Asterisks are numerical 
results, blue and red curves correspond to curves in upper 
panels. 



tions: 



-e 2 [|ai| 2 sin 2a; + |a 2 | 2 sin(2x - tt/2)] 



ai = iAidi + 1 + iai cos x — ka^i 
a 2 = «A 2 a 2 + R + ia 2 cos [x — 7r/4) 



KAO-2- 



(2) 



The dynamics for a given R < 1 depends only on the 
the scaled driving e 2 = C,E\ where C, = ^gs, two scaled 
detunings Ai. 2 = 8^2/A and a scaled damping ka = 



^ K ' A ' ; all scaled frequencies (including cooling rates) are 
given below as a fraction of A. 

The experimentally adjustable parameters are e, both 
the detunings Ai t 2 and R. We assume ka — 1, though 
the analytical expressions are for arbitrary ka- Varying 
driving power V ~ 1 — 10 mW , but leaving the cav- 
ity/nanosphere properties unchanged means A remains 
constant, but e 2 varies from ~ 1 — 100. 

Following the usual procedure, we replace operators by 
their expectation values and linearise about equilibrium 
fields by performing the shifts a\ — > ai+ai, 02 — > 0:2 + 0.2 
and x —¥ xo + x. Hence we find equilibrium photon fields, 
cei = [ka — i&iY and ot2 = R [ka — iAf ] _ as well as 
position tan 2:ro = |a2| 2 /|ai| 2 . 

Here, A? = Ai + |(1 + cos2x ) and Af = Ai + |(l + 
sin2xo). The dimensionless mechanical frequency is: 



J M 



(Ai,A 2 



= 2e 2 (|ai| 2 cos2a;o 



|a 2 | 2 sin2x ). (3) 



Closely related forms of this "self-trapping" frequency 
expression have been noted previously |14l - [l7J but the 
implications, other than for xq — 0, have not been inves- 
tigated. 

To first order, the linearised equations of motion are: 

x = — ujf.jx — e 2 (<7i sin2xo — (72 cos2xo) 
di = i A1CL1 ~ iaixsui2xo — kacli 
a.2 = iA^a?, + ia^x cos2xq — KA&2 (4) 

where gi — (a*Oi + ona*). From the above, we can ob- 
tain the contribution from the two photon fields to the 
optomechanical cooling: 



r 

2 

where 



e 2 KA 

2lom 



[S\{um] + S 2 {u:m) - Si(-Wm) - &(-Wm)] 

(5) 



|ai| 2 sin 2 2x 
Sl(w)= [Af- W ]» + «^ ;S2(w) 



|a2| 2 cos 2 2a;o 



(6) 



(net cooling occurs for T < 0). We also calculate a nu- 
merical r by evolving the equations of motion in time 
and looking at the decay in x(t) (its variance in partic- 
ular). The analytical cooling rates give excellent agree- 
ment with numerics in all but the strongest cooling re- 
gions. 

The single-field cooling resonance r2 occurs for xq ~ 0, 
Af ~ and Af = -Wm, thus here for Ai = —1 (i.e. 
8\ = —A). Conversely, there is also a single-field cool- 
ing resonance rl for xq — vr/4 (note that since we con- 
sider only the case R < 1, i.e. field 2 is always driven 
more weakly than field 1, the latter situation does not 
correspond simply to an interchange in the role of the 
fields). Away from these extreme cases, both cooling res- 
onances are split by the effect of the other field, wherever 

< Xq < 71-/4. 

From Fig 12] we see r2± occur for the same A2, thus 
the same equilibrium photon field 02; however they cor- 
respond to photon fields a\ and a\ respectively and thus 



to different A^ = ±yi, where 2y\ is the splitting (about 
A^ = 0) between r2+ and r2— seen in Fig(2] 

The transformation ai — > a\ leaves both the mechan- 
ical frequency ui m (Af, A2) and cco(Aj = ,A2) unchanged. 
Hence the cooling rates are similar for both r2±. 

We can estimate the splitting A* by requiring 



UM 



(A+,A 2 )- WM (Ar,A 2 )~A? 



(7) 



since wmCA^, A2) — A 2 are the conditions for the op- 
tomechanical resonance r2±. From Eqs. flTB)) and © we 



see: 



±2/1 = ±1 



2c 2 



(Af) 2 cos 2x 



4- V2e 
^ A 2 +A/2- 



(8) 



Similarly, j/2, 



2e 2 R 2 _ K 2 

(AJ) a sin2£c A- 



Close to r2, we can simplify y\ 

the splitting between rl± is ±2/2 = ± 

Thus the splittings increase with driving power and R. 
While in Fig[2j corresponding to R — 0.5, three reso- 
nances (r2— , rl— and rl+) overlap, in Fig[3j for R = 1 
rl± are well separated and the double resonance involves 
only r2— , rl — . 



Cooling Rate R = 1 




FIG. 3: (Colour online) Cooling rates as a function of scaled 
detunings. For R = 1 the behaviour of mode 1 and mode 
2 is equivalent, thus a high degree of symmetry is evi- 
dent. The splitting of the resonances is much larger, but 
very strong cooling maximum (dark blue) at the double- 
resonance of r2— , rl— is seen, with a matching heating max- 
imum (white/orange) for Af 2 > 0. 

We now analyse the relative merits of single-resonance 
versus double-resonance cooling. Single-field cooling cor- 
responds to r2 in FigJ5] Cooling rates are obtained from 
Eq|5]by taking x ~ 0, A* = and Af = -w M - This 
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FIG. 4: Comparison between single-resonant (r2) and dou- 
ble resonant (DR) i.e. simultaneous rl— and r2— resonant 
cooling rates as a function of laser driving power (R = 0.5), 
showing that the double resonance corresponds to two orders 
of magnitude greater cooling for strong driving. Inset shows 
cooling maximum of r2 at weak driving e 2 ~ 1/2; here the sys- 
tem is on the edge of the sideband-resolved regime ujm = ka- 
In contrast for DR, at P — WmW , ujm/ka — 4 and the rl — 
and r2— sidebands are very well resolved. 



regime was investigated in [15jj and Eqj5] reduces to ex- 
pressions therein (in unsealed units). However, we can 
give good approximations to cooling rates purely in terms 
of experimental parameters (driving power, R and k). 
Assuming S2(-wa;) S> S^+^a/) and that the field 1 
contribution to cooling is negligible, near r2, the me- 
chanical frequency ui M = ^-. Hence, as shown in the 

appendix, the Single Resonance (SR) cooling rate be- 
comes: 



SR 



V2 



-(2 £ 2 + 4r 



(9) 



(recall this is a scaled cooling rate thus given in units of 

Single-field cooling is a maximum if e = k? a J \/2 where 
ljm = k/2 (in unsealed units) and is thus at the edge of 

the resolved sideband regime. Here, — Tsr ~ *%■', this 
gives optimal cooling T ~ k only if R ~ 1. This cooling 
maximum is independent of ka'- it depends only on R. 
As the driving is increased, if 2e 2 3> k a , 



SR 



R 2 k\ 

2^/2t 



oc 1/e 



(10) 



rapidly as the driving amplitude is increased: the cooling 
cannot be improved by increasing the driving amplitude. 

To obtain the corresponding double-resonant rate 
(DR) one must first identify the e-dependent pair of de- 
tunings for which — wm(Ai,A2) ~ Af w A|. Even if 
rl — ,r2— do not cross in FigJ5J the sidebands approach 
within their width k and overlap significantly. One can 
still obtain a good approximation to the cooling rate in 
terms of driving parameters (see Appendix) . In this case, 
there are contributions to cooling from both field 1 and 
field 2. Adding them both, 



DR 



e 2 {R 2 +R 4 



(kAUm(Wm + K, A ) 



(11) 



where the frequency is given by the expression, ui 2 M = 



— Sp- + i \/k a + 8e 2 . The contribution from mode 2 ex R 2 
while that of mode 1 ex i? 4 ; both contribute significantly 
for R > 0.5. Assuming wm 3> KA, this reduces to 
-T DR « 2- 3 / 4 (i? 2 + i? 4 )e 1 /2/ KA , hence T DR ex e 1 / 2 . 

FigH] shows that Eqs|9]and Eqs[TT]both give excellent 
agreement with exact numerics. In the double-resonant 
case, the cooling is stronger and increases with increasing 
e. In contrast, the single- resonant cooling Tsr oc 1/e and 
this cannot be improved by increasing e. Self-trapping 
cooling cannot be considered simply in terms of an addi- 
tive contribution from two intracavity intensities in Eqj6] 
the response of u>m to the driving is also important. In 
the double-resonant case, ujm oc e 1 ' 2 for strong driv- 
ing. In contrast, for the single-resonant case u>m oc e 
and strong driving pushes the r2 resonance into the far- 
detuned regime. A study of the quantum cooling shows 
that the minimum phonon numbers attainable n m in "C 1 
for single and double resonant cooling for strong driving 
e 2 ^> 1. For weak driving, strong cooling can be obtained 
with single field cooling, but at the edge of the resolved 
sideband regime uim = k/2, less favourable for ground 
state cooling. 

Conclusion Our study shows that the two-mode self- 
trapping regime has a raft of of unexpected features, 
including the split side-bands, strong cooling at blue- 
detuning and simultaneous heating and cooling reso- 
nances. Although other proposals also permit strong 
cooling rates, the multiple sidebands provide an excep- 
tionally broad region of strong cooling, offering consid- 
erable robustness to experimental errors in the driving 
power, detunings and even the phase (a variation of </> of 
order 30% will not appreciably perturb the strong cool- 



Thus the single resonance cooling rate falls off quite ing) 
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I. APPENDIX 

A. Minimum phonon number 

From quantum perturbation theory we can show that 
Rn^m, the rate of transition from state n to n+1 is: 



Rn^n+l = (n + 1) 



Similarly: 



£ 2 ka 



(ftW + S 2 (u M )) (12) 



-Rrwn-l = U (Sl(-WAf) + Sa(-WM)) (13) 

For n » 1, this gives the cooling rate of (EqJS]). How- 
ever, with the exact expressions we can also show the 
equilibrium mean phonon number to be: 



Si(-ujm) + S 2 {-u) M ) 

mm Si(um) + S 2 {ujm) - Si(-lum) - S 2 (-ojm) 

(14) 
The parameters for strongest cooling do not necessar- 
ily yield the minimum phonon numbers, because the rate 
depends upon the difference between the optical heating 
and cooling, while the phonon number depends upon the 
ratio between the two. In Figs|5]and[6]we present colour 
maps comparing the cooling and minimum phonon num- 
bers for both R — 0.5 and R = 1 respectively. Fig(5] cor- 
responds to the same parameters as FigJU It shows the 
broad strong cooling region corresponding to the three 
distinct cooling resonances r2— , rl± which are only par- 
tially resolved. The result is a strong cooling region of 
about 1MHz width, providing the advantage of a strong 



cooling regime r ~ k insensitive to experimental detun- 
ings. For Fig|5]in contrast, the map has a high degree of 
symmetry, since the role of the modes is interchangeable; 
the larger R means the splitting is larger, so rl± are fully 
resolved at the double resonance: thus only r2— , re- 
contribute to the maximal cooling region. Nevertheless, 
this is the point where the strongest cooling (r = 0.95K.4) 
is obtained. 
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FIG. 5: Colour online) maps of cooling rate and minimum 
phonon number for parameters R — 0.5 equivalent to Fig[2] 
The cooling rate map shows information very similar to the 
shifted curves in Fig[2] The white lines indicate the locus of 
the single field resonances rl± (where Af ± = ujm and r2±, 
where (Af± = ujm)- Even at the Double Resonance (where 
the two white lines intersect), there is still a contribution from 
the partly resolved third sideband of rl+, giving a very broad 
strong-cooling region (dark blue). The corresponding broad 
strong-heating region is also seen for positive AJ. NB: The 
axis correspond to corrected detuning scaled to A = re, Ai j2 , 
not the experimental <5i,2 of Fig[2] R = 0.5 maximum cooling 
is T = -0.59. 



The figures show that the largest cooling rates are 
found in the symmetric double resonance regime (the 
double resonance for R = 1), the phonon occupancy 
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B. r2 and single field cooling 

The single-field cooling rates are obtained from Eq[5] 
by taking x ~ 0, A^ = and Af = —u>m- 
Hence, 



r/2 



e 2 ^ , |2 
ad 



2w 



M 



1 



1 



\U>M 



^21 



Wm 



Af] 2 



We can obtain the precise form given previously 
J15J if we replace |a 2 | 2 using EqUU |^4j ~ 2x ~ 

.2 

7-7 — , j-t|y2 i — a- and revert to unsealed units. 

The above, in principle, requires a full numerical solu- 
tion of the equations: 



fli 
tan 2xq 



[KA-iAf] 1 ; a 2 ^R[K A -iA^Y 



l« 2 | 2 /l«il 2 , 



(16) 



to find equilibrium positions Xq and fields. However 
we can obtain a good approximation in closed form us- 
ing experimental parameters (driving power, detunings). 
Hence, 



2e 2 



J M 



(17) 



For resonant cooling S%(— Um) ^ S2(+wm) and cool- 
ing/heating by field 1 is negligible. So, 



r 



e 2 |a 2 | 2 
wmM 



(18) 



FIG. 6: As for Fig[5]but for R — l. The fields are equivalent, 
thus the cooling/heating maps have a high degree of reflection 
symmetry. The white lines indicate single field resonances, 
double resonance occurs where the two white lines intersect. 
In this case, the splitting between rl± is larger so the double- 
resonance involves only rl— and r2— . Fhe R — 1.0 max 
cooling is F = -0.95, Af = Af = -2.72. 



is slightly lower in the usually analysed [15| single-field 
regime where one strong field, at near zero detuning, 
traps while the other field cools. Nevertheless, the mean 
equilibrium occupancy is very small in both cases (less 
than a tenth of a phonon). The above analysis of mini- 
mum phonon occupancy has not considered other sources 
of heating (photon scattering, background gas collisions). 
The effect of these other heating rates has been analysed 
in [15[ where ground state cooling is shown to be achiev- 
able for optimal cooling rates (r ~ «). Thus the 1-2 order 
of magnitude increase in cooling in the double-resonance 
region means that this is the most favourable regime. 



and |«2 1 2 = 



—5-, hence the Single Resonance (SR) 



cooling rate becomes: 



SR 



R 2 en 2 A 



V2(2e 2 + k\) 



(19) 



NB: this is a scaled cooling rate thus given in units of A. 
This gives a maximum cooling rate if e = k a /\/2 where 



SR 



4 



(20) 



It is worth noting that this maximum is independent 
of ka'- it depends only on R. (Within the underlying 
assumption that u>m » ka) Even for R = 1/2 , the 
maximum cooling is Tsr/2 ~ 1/32 far from optimal 
r o pt/2 ~ 1. 

As the driving is increased, if 2e 2 ^> k\, 



V SR 



R 2 k\ 

2\/2e 



ex 1/e 



(21) 



Thus the single resonance cooling rate falls off quite 
rapidly as the driving amplitude is increased: the opti- 
mal cooling cannot be attained by increasing the driving 
amplitude. 



C. r2 and rl — overlap: double-resonance cooling 

As illustrated in Fig|2] the resonances r2— and rl— 
never actually coincide: they undergo something remi- 
niscent of an avoided crossing (recall that Figgis a map 
of the classical cooling, not of underlying quantum eigen- 
values. Nonetheless, the similarity between the classical 
linearisation and the effective Hamiltonian for the quan- 
tum fluctuations has a very similar structure. 

The double-field cooling rates for the region of closest 
approach are estimated from Eq[5] by assuming that the 
two resonances actually cross, in other words both fields 
are, to a good approximation, simultaneously resonant. 

The first task is to identify the e-dependent pair of 
detunings for which: 



w M (Ai,A 2 )«A? 



AX 

^2 



(22) 



One can search numerically for detunings which give 
near-simultaneous resonances. However, for moderate 
sin2xo — tan 2xo — 2xo, a closed form can be obtained. 
Provided R is small the double resonance still falls in the 
small angle regime (including our FigJ2J. In this case, 
From Eqll6l 



2x n ~R 



2 (AT) 2 + 4 
(A|) 2 + k\ 



R 1 



(23) 



so if AT = A2 = — u>m we obtain xq ~ R 2 /2. 
The mechanical frequency: 



lj 2 m = 2e 2 \a 1 \ 2 



2c 2 



2e 2 



(AT) 2 
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Eg s l26l and [27l should be contrasted with the behaviour 
of the singly resonant cooling EqllOl In the double- 
resonant case, the cooling increases, without limit, as 
a function of e. For R = 1/2 and e 2 ~ 100 (ie 8mW 
power), r ~ 0.4 A In Fig|4] the approximate expres- 
sions Eqs ll0|26|27l are shown to give quite good agree- 
ment with cooling rates obtained from numerical solution 
of the equations of motion, without linearisation. 

The very large damping rates provided by the double- 
resonant regime have the added advantage of relative in- 
sensitivity to the initial preparation. Although the anal- 
ysis assumes small oscillations about equilibrium xq, a 
resonator prepared at the antinode of the trapping field 
(x = 0) is very rapidly pulled towards equilibrium and 
oscillates about x = xq. Consequently there is no need to 
provide any initial displacement. In addition, the multi- 
ple resonance region shown in Figgis surprisingly robust 
to errors in the relative phase between the two modes. In 
FigJS]we see that a phase error of order 20% makes rel- 
atively little difference (and may even enhance the cool- 
ing). 



D. Symmetric double-resonance cooling 

The best regime for strong cooling is one in which the 
fields are equally strong and where both are resonantly, 
red detuned. In this region: 
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J M 
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(28) 



this means: 



J M 



(25) 



In this case, there are contributions to cooling from 
both field 1 and field 2. Adding them both, 



and: 
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In the limit lom » ka' 



Mf 



(29) 



DR 



e 2 (R 2 +R i ) 

«AWm(Wj, + K 2 A ) 



(26) 



J- opt — 



2" 9 / 8 Vi 



U) M 



(30) 



where we will substitute Eq l25l to evaluate the frequency 
ujm. As ever, u M > ka- 



r_Di?. 
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This is not simply a special case of Eqf27] because we 
have moved away from small xq into the regime x = n/8. 
However the difference between Eq[27] and the above is 
less than a factor of 2. 




-4e+05 



FIG. 7: Shows that the strong cooling region where the split 
sidebands r2— ,rl— ,rl+ overlap is reasonably insensitive to 
the phase difference between the two modes. A variation of 
about 30% about the usual dephasing of it/2 is tolerable (and 
can even give enhanced cooling). 



